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Within the low-energy effective Minimal Supersymmetric extension of Standard 
Model (effMSSM) we calculated the neutralino relic density taking into account 
slepton-neutralino and neutralino-chargino/neutralino coannihilation channels. We 
performed comparative study of these channels and obtained that both of them 
give sizable contributions to the reduction of the relic density. Due to these coan- 
nihilation processes some models (mostly with large neutralino masses) enter into 
the cosmologically interesting region for relic density, but other models leave this 
region. Nevertheless, in general, the predictions for direct and indirect dark mat- 
ter detection rates are not strongly affected by these coannihilation channels in the 
effMSSM. 
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I. INTRODUCTION 



Measurements of the cosmic microwave background radiation |L| imply that the universe is 
almost flat fl = p/p c — 1, where p c — 3H 2 /8itGn is the critical closure density of the universe, 
Gm is Newton's constant and H = 100 h km/sec/Mpc is the Hubble constant with h = 0.7 ± 0.1 
£f) ' ||. A variety of data ranging from galactic rotation curves to large scale structure formation 

and the cosmic microwave background radiation imply a significant density 0.1 < flh 2 < 0.3 || 
of so-called cold dark matter (CDM). It is generally believed that most of the CDM is made 
of weakly-interacting massive particles (WIMPs) j|. A commonly considered candidate for the 
WIMP is the lightest neutralino, provided it is the lightest supersymmetric particle (LSP) || in the 
Minimal Supersymmetric extension of Standard Model (MSSM). Four neutralinos in the MSSM 
being mass eigenstates are mixtures of the bino B, wino W and higgsinos H®, H®, and the LSP 
can be written as a composition x = Xi = N\\B + Ni 2 W + Ni 3 H® + AT 14 iJ°, where Nij are the 
entries of the neutralino mixing matrix. In SUSY phenomenology one usually classifies neutralinos 
as gaugino-like (with V ps 1), higgsino-like (with V ~ 0), and mixed, where (gaugino) purity is 
defined as V = |A n | 2 + |7V 12 1 2 . 

In most approaches the LSP is stable due to R-parity conservation |6|. The neutralino, being 
massive, neutral and stable, often provides a sizeable contribution to the relic density. The contri- 
bution of neutralinos to the relic density is strongly model-dependent and varies by several orders of 
magnitude over the whole allowed parameter space of the MSSM. The neutralino relic density then 
can impose stringent constraints on the parameters of the MSSM and the SUSY particle spectrum, 
and may have important consequences both for studies of SUSY at colliders and in astroparticle 
experiments [Q. In light of this and taking into account the continuing improvements in deter- 
mining the abundance of CDM, and other components of the Universe, which have now reached 
an unprecedented precision one needs to be able to perform an accurate enough computation 
of the WIMP relic abundance, which would allow for a reliable comparison between theory and 
observation. 

In the early universe neutralinos existed in thermal equilibrium with the cosmic thermal plasma. 
As the universe expanded and cooled, the thermal energy is no longer sufficient to produce neu- 
tralinos at an appreciable rate, they decouple and their number density scales with co-moving 
volume. The sparticles significantly heavier than the LSP decouple at the earlier time and decay 
into LSPs before the LSPs decouple themselves. Nevertheless there may exist some other next-to- 
lightest sparticles (NLSPs) which are not much heavier than the stable LSP. The number densities 
of the NLSPs have only slight Boltzmann suppressions with respect to the LSP number density 
when the LSP freezes out of chemical equilibrium with the thermal bath. Therefore they may 
still be present in the thermal plasma, and NLSP-LSP and NLSP-NLSP interactions hold LSP in 
thermal equilibrium resulting with significant reduction of the LSP number density and leading to 
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acceptable values even with a rather heavy sparticle spectrum M. These coannihilation processes 
can be particularly important when the LSP-LSP annihilation rate itself is suppressed. 
The number density is governed by the Boltzmann equation |j 10 



^+3Hn = -(av)(n"-nl q ) (1) 

with n either being the LSP number density if there are no other coannihilating sparticles, or 
the sum over the number densities of all coannihilation partners. The index "eq" denotes the 
corresponding equilibrium value. To solve the Boltzmann equation one needs to evaluate the 
thermally averaged neutralino annihilation cross section {<jv). Without coannihilation processes 
(av) is given as the thermal average of the LSP annihilation cross-section a xx times relative velocity 
v of the annihilating LSPs 

(av) = {<T xx v), (2) 

otherwise, it is determined as (av) — (a c sv), where the effective thermally averaged cross-section 
is obtained by summation over coannihilating particles |^ 10 



oq eq 

. . ^— \. \ n i n j 

y 

If no denotes the nowadays number density of the relics, the relic density is given by 

n = ^. (4) 

Pc 

Many increasingly sophisticated calculations of the relic density of neutralinos in supersym- 
metric models, with various approximations both in the evaluation of (a c ffv) and in solving the 
Boltzmann equation were performed . Following §|] we 

briefly remind below of the major developments in the field. 

Perhaps for the first time strong constraints for pure photino relic abundance were obtained 
in Jl2]| . The first analysis of the general neutralino case was performed in Several other 

early papers subsequently appeared with more detailed and elaborate analyses. The two-neutralino 
annihilation into the ordinary fermion-antifermion (//) final states through the Z-exchange was 
computed in detail Jig ], The first complete analysis of the neutralino annihilation into W + W~ , 
Z°Z° and Higgs-pair final states was performed in |lq| . The Higgs contribution into neutralino 
annihilation was first computed in [ ^6|p^ ] . For the pure gaugino-like and higgsino-like neutralinos 
(where several important resonances and final states are absent) all the annihilation channels 
were considered in ]jfi| |. A first complete set of expressions for the product of the cross section 
times velocity using the helicity amplitude technique was computed in fisfl . The relic density 
calculation was made by expanding the (av) as a power series in neutralino velocity. The angular 
and energy integrals in such a case can be evaluated analytically and the remaining integration 
over temperature was performed numerically. When expanded in the nonrelativistic limit, these 
give expressions for the first two coefficients of the partial wave expansion. In the early papers the 
partial wave expansion of the (av) was used in most cases. The method is normally expected to 
give an accurate enough approximation in many regions of the model parameter space because the 
relic neutralino velocity is expected to be highly non-relativistic. However, it fails near s-channel 
resonances (quite high-energetic ) and thresholds for new final states, as was first pointed out in p0[ | 
and further emphasized in P,ETU23]. In particular, it was shown [plf that due to the very narrow 
width of the lightest supersymmetric Higgs h, in the vicinity of its s-channel exchange the error can 
be as large as a few orders of magnitude. Therefore a relativistic treatment of thermal averaging 
is necessary. A recent detailed analysis (3(| showed that in the case of the often wide s-channel 
resonance exchange of the pseudoscalar Higgs boson A, the expansion produces a significant error 
fl9[l . Furthermore subdominant channels and often neglected interference terms can also sometimes 
play a sizeable role. 

The proper formalism for relativistic thermal averaging was developed in || , and used in f2~if| . 
A more accurate treatment of the heat bath for both annihilating particles involving two separate 
thermal distributions was considered in flp3fl. A very useful compact expression for (av) as a 
single integral over the cross section was for the first time derived in M. The DarkSusy code 
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where the relic density of neutralinos is numerically computed without the partial wave expansion 
approximation was developed in [fj"l| . 

An additional very strong reduction of the relic abundance of WIMPs through coannihilation 
was first discovered in pcj. There are regions in the MSSM parameter space where higgsino-like 
LSP, light chargino and next-to-lightest neutralino masses become nearly degenerate and all three 
species can exist in thermal equilibrium. Their mutual coannihilation is often important, and 
even dominant j^j3j],[l7J] . In the coannihilation with the LSP can be involved any SUSY particle, 
provided its mass is almost degenerate with the mass of the LSP In the low -energy effective 

MSSM (effMSSM), where one ignores restriction from unification assumptions and investigates the 
MSSM parameter space at the weak scale [^0| , p6|j4C[| there is, in principle, no preference for the 
next-to-lightest SUSY particle. Nevertheless due to quite reasonable and commonly used sets of 
free parameters, when all gaugino mass eigenvalues are calculated in terms of entries of gaugino 
mass matrices (fi, M2, tan/3), the coannihilations between gauginos are expected to occur most 
often, in the effMSSM ||,||[l0|. 

The relativistic thermal averaging formalism [^| was extended to include coannihilation processes 
in JlfH , and was implemented in the DarkSusy code [Til ] f° r coannihilation of charginos and heav- 
ier neutralinos. In was found |l(J that for higgsino-like LSP such a coannihilation significantly 
decreases the relic density and rules out these LSPs from the region of cosmological interest. For 
the highly bino-like LSPs, the reduction of the relic density due to the coannihilation is not strong 
enough to avoid an overdosing of the universe. 

The importance of the neutralino coannihilation with sferminos was emphasized and investigated 
both for sleptons PJ29|] or stops p3| , p4| and sbottoms | 3lfl in the so-called constrained MSSM 
(cMSSM) f^^fTjor in supergravity (mSUGRA) models~|rj] . 

In the most popular mSUGRA model SUSY breaking occurs in a hidden sector and is 
communicated to observable sectors via gravitational interactions. The model has a minimal set 
of parameters: mo, mi/2, Ao, tan/3 and sign(yu). Here mo is the universal scalar mass, is the 
universal gaugino mass and Aq is the universal trilincar mass, all evaluated at Mqut, tan/3 is the 
ratio of Higgs field vacuum expectation values and /1 is a Higgs parameter of the superpotential. In 
particular, within the framework of the mSUGRA, it was found ^4|,^8| that at large tan /3, indeed 
large new regions of model parameter space gave rise to reasonable values for the CDM relic 
density, due to off-resonance neutralino annihilation through the broad A and H Higgs resonances. 
There are strong correlations of sfermion, Higgs boson and gaugino masses in mSUGRA originating 
from unification assumptions. In regions of mSUGRA parameter space where x an d fi (or other 
sleptons) were nearly degenerate (at low mo), coannihilations could give rise to reasonable values of 
the relic density even at very large values of m 1( / 2 , at both low and high tan/3 |p|,pT|. In addition, 
for large values of the parameter Aq or for non-universal scalar masses, top or bottom squark 
masses could become nearly degenerate with the %, so that squark coannihilation processes can 
become important as well [[53 34 1. Therefore due to slepton and squark coannihilation effects, 



the relic density can reach the cosmologically interesting range of 0.1 < Vlti 1 < 0.3 jj). In the 
mSUGRA LSP is naturally almost pure bino-like as was first noticed in from the point of view 
of low-energy SUSY and CDM. 

Having in mind investigation of future prospects for direct and indirect detection of LSPs we 
follow the most phenomenological (general) view, not bounded by theoretical restrictions from 
sfcrmion/gaugino/Higgs mass unifications, etc. To this end we need maximally general and ac- 
curate calculations of the relic density within the low-energy effective MSSM scheme (effMSSM) 
p^ , [fO[ | . The only available high-level tools for these calculations was the DarkSusy code (the best 
code to our knowledge at the moment, this paper was started). Unfortunately the code calculates 
only neutralino with neutralino/chargino coannihilations (NCC), which is not sufficient in the case 
of bino-like LSPs, when neutralino-slepton coannihilation (SLC) and neutralino-squark coannihi- 
lation are claimed to be dominant |^,^9 33 3^ , |3l|| . This paper is aimed at a comparative study of 
NCC and SLC channels, exploration of relevant changes in the relic density and investigation of 
their consequences for detection of cold dark matter particles in the effMSSM. 



II. THE effMSSM APPROACH 

As free parameters in the effMSSM we use the gaugino mass parameters Mi, M2; the entries 
to the squark and slepton mixing matrices n^,m%,m%,m 2 -,m 2 ~ for the 1st and 2nd generations 

U D l\ L 
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and m Q 3 > m f, m %i m \ 3 > rn \ 3 f° r the 3rd generation, respectively; the 3rd generation trilinear soft 
couplings A t ,Ab,A T - the mass itia of the pseudoscalar Higgs boson, the Higgs superpotential 
parameter /x, and tan/3. To reasonably reduce the parameter space we assumed m 2 ^ = m 2 ^ = m 2 ^', 

m% = m 2 ~ = m% ; m 2 ~ = rn 2 - ; m 2 - = m\ and have fixed Aj, = A? = 0. The remaining parameters 

1 B Q3 H L H3 L3 

defined our effMSSM parameter space and were scanned randomly within the following intervals: 

-1 TcV < Mi < 1 TeV, -2 TcV < M 2 , fi,A t <2 TeV, 1.5 < tan^ < 50, 
50 GeV < M A < 1000 GeV, 10 GeV 2 < m|, m\, m| 3 , m| 3 < 10 6 GeV 2 . 

We have included the current experimental upper limits on sparticle masses as given by the Particle 
Data Group [Q. As in ficfl , we have used the limit nk± > 85 GeV for the lighter chargino mass. 
The current limits on the rare b — > 57 decay [ ]45[ and the g — 2 limits |4q ] have also been imposed. In 
agreement with a flat accelerating universe ||]", we assume 0.1 < Vthr < 0.3 for the cosmologically 
interesting region. The calculations of the ncutralino-nucleon cross sections, and direct and indirect 
detection rates follow the description given in |5]j4C||. 

We have evaluated the relic density of the LSP under ignoring of any coannihilation (IGC), 
taking into account only NCC or SLC separately, as well as including both coannihilation channels 
(BCC). To this end in our former code |l0| DarkSusy procedures of (a c gv) evaluation and solution 
of Boltzmann equation were implemented. All coannihilations with two-body final states that 
occur between neutralinos, charginos and sleptons, as long as their masses are m, < 2m x were 
included. The Feynman amplitudes for NCC and SLC were taken from DarkSusy |ll| and |47|]|], 
respectively. We calculated {(J c sv) and flh 2 following the relevant DarkSusy routines H], which 



wc have merged with code 47 in a way that guarantees the correct inclusion of SLC. 

If all sleptons, neutralinos and charginos in question are substantially heavier than the LSP 
(rrii > 2m x ) and no way for coannihilations, the resulting relic density flh 2 = VlK^qq = VlK^qq — 
r2ft,| LC is equal to Q x h 2 of former results obtained without any coannihilations (with (a e f{v) = 
(a xx v)). When at least one of coannihilation channels (NCC or SLC) is indeed relevant, the 
^igc (ignorance of any coannihilation) is calculated with 

f n 7\ 2 

(o-eff^)iGC = (crxx v ) [ I ' ( 5 ) 



where n eq includes all open coannihilation channels. This formula allows a comparative study 
of results, relevant to one or both coannihilation channels, always delivering a decreasing ratio 
Qhg OA /nh 2 GC < 1 in accordance with |loj and sometimes contrary to |To|,p|]. We introduced 
^^coa as a common notation for f2/ig CC , f2/ij^ cc or flh^ LC . 



III. RESULTS AND DISCUSSIONS 

We performed our calculations in the effMSSM approach given above and results of our consid- 
erations (scatter plots) are presented in Figs. |lj — |^ and Figs. |7]-|| for neutralino relic density and 
CDM observables, respectively. 



A. Relic density 

The general view of the reduction effect on the relic density (RD) due to NCC, SLC and BCC 
are shown in Fig. [l] as ratios f2/i 2 ; , OA /r2/i 2 GC together with comparison of NCC against SLC in the 
form of the ratio ^^|lc / ^^-ncc ■ ^ n the- basis of our sampling (50000 models tested) the maximum 
factor of RD decrease due to NCC is about 2 • 10 -3 for m x w 200 GeV, while SLC reduces the 
RD maximally by a factor of 8 • 10~ 4 for m x w 300 GeV. Both reduction factors are roughly of 
the same order of about 10 -3 . These results depend on the applied experimental limits on the 
second-lightest neutralino, chargino and stau masses. If there were no limits implemented on their 
masses, the factor of relative RD reduction due to NCC could reach a maximum value of 2- 10 -5 for 
models with m x « 40 GeV. But in our case, the current experimental limits are m x ± > 85 GeV 
and rrif > 81 GeV, and therefore the critical LSP mass that enables non-negligible NCC and SLC 
contributions is also of the same order (m x > 80 GeV). 
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FIG. 1. Effects of neutralino-chargino/neutralino (NCC) and slepton-neutralino (SLC) co-annihilations. 
Panels a)-d) display ratios S7/j,b CC /fi/if GC , £lh^ cc /flhi GC , flhg hC /{lh 2 GC , and flh^ cc /£lhg hC , respectively. 
The maximal reduction factors for both NCC and SLC are of the order of 10~ 3 . 




ffl 



700 800 
mLSP, GeV 



FIG. 2. The same as in Fig. [l] a), b) and c), but plotted together. Here £lh 2 iGG /Qh 2 an, Qh 2 



[(;<'■ ^"NCc/^IGO 



and flhg LC /Q.h 2 GC are marked with squares, triangles, and stars, respectively. Therefore, a square filled 
with a star (triangle) depicts a model that is only affected by SLC (NCC), while the other coannihilation 
channel in the majority of models gives negligible contribution. 
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From Fig. || one can see that the reduction of RD by coannihilations is mainly due to either 
NCC or SLC. The other channel of coannihilation plays no role or leads only to a much smaller 
further reduction. Although other coannihilation processes besides NLSP-LSP can in principal 
occur (including the next-to-NLSP (NNLSP) and next-to-NNLSP, etc), Fig. | demonstrates that 
a stau f as a NLSP indeed entails a dominant SLC effect, while a next neutralino \2 or chargino 
X ± as a NLSP indeed entails a dominant NCC effect. 




A A 



FIG 



3. Ratio fi/iNCc/^SLC versus m x 
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Stars indicate that the f is the NLSP, triangles mean that the 



light chargino \ is the NLSP, small filled squares mark the models where the second-lightest neutralino 



\2 is the NLSP. One sees that if f is the NLSP, the SLC necessarily dominates, while \2 or x 
NLSP always leads to dominant NCC. 
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From Fig. [| one can notice that mass differences rrif — m x < 20 GeV lead to a RD reduction factor 
of 0.5-0.005; 20 GeV < m f - m x < 40 GeV lead to the factor of 0.8-0.01 and mass differences 
40 GeV < mf — m x < 100 GeV can still lead to factors smaller than 0.3 due to SLC. Mass 
differences m X2X ± — to x < 2 GeV lead to a RD decrease by factors of 0.1-0.005; 2 GeV < m X2 ^ x ± — 
m x < 30 GeV lead in general to factors 0.9-0.02 due to NCC. For both kinds of NLSPs, the 
coannihilation effect may become negligible if tonlsp — % > 30 GeV, and necessarily becomes 
negligible if tonlsp — to x > 100 GeV. Therefore, future increase of the lower mass limits for all 
possible NLSP (at Tevatron or LHC) can, in principle, strongly reduce the importance of the effect 
of any of the coannihilation channels. 

Although we have implemented the coannihilation opening threshold of rrii — 2m x , it was 
found that a SLC-reducing factor less than 0.5 (0.1) occurs only for mf < 1.12 (1.05) m x . Ac- 
cordingly, a NCC-reducing factor less than 0.5 (0.1) appears for m x ± < 1.16 (1.10) m x and 
m X2 < 1.11 (1.08) m x . Therefore for all channels of coannihilation, relevant effects occur if the 
mass difference between the coannihilation partner and the LSP is within 10-15%. This is in an 
agreement with previous considerations |2^j8] ^ , |l^ , ^p9lj3^j3l|l . 

From Fig. || one can also see that charginos and neutralinos come to lie close in mass to the LSP 
more often than staus (and other sleptons). This is an expected result of correlations in the gaugino 
sector of the effMSSM, as mentioned in the Introduction, which explains the NCC dominance over 
SLC seen from Figs. || and |[ If one manages to construct a SUSY model where the LSP mass 
is almost always degenerate with one of the slepton masses (as for example, in mSUGRA models 
with bino-like neutralinos) the dominant coannihilation channel will be SLC. 
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FIG. 5. Illustration of the shifting of effMSSM models inside and outside the cosmologically interesting 
range 0.1 < Q.h COA < 0-3 due to NCC and SLC. Models with Qhi GC , f2ftg LC , Qh^cc an d ^bcc are marked 
with empty circles, filled circles, small dots, and empty squares, respectively. Therefore, a superposition 
of all those symbols corresponds to a model which is totally untouched by coannihilation. A black-framed 
filled circle marks a model which is untouched by SLC (£l/i| LC = Q,hf GC ), but shifted down due to NCC. 
If the corresponding f2/i| cc (which is equal to Qh^ cc ) remains within this range, it still presents in the 
figure below the black-framed filled circle as an empty square with a black dot inside. By analogy, an 
empty square with a filled circle inside gives a model which was shifted into the region due to SLC only 
(f2/i| cc = Q.hg LC ), and if the corresponding Qhi GC also is in the cosmologically viable range, it is located 
above the filled square as an empty circle with a dot inside. One can notice that a quite big amount of 
models is shifted out of 0.1 < fi/i 2 < 0.3 due to NCC (grey circles). 

In Fig. |^ all calculated relic densities (f2/ij GC , f2/i§ LC , f2/i^j CC and f2/i| cc ) are depicted in the 
cosmologically interesting region 0.1 < f2/i COA < 0.3. There is a quite big amount of models 
(mostly with lower LSP masses) which are completely unaffected by coannihilation. When at least 
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one of coannihilation channels is relevant, the RD decreases and some cosmologically unviable 



models with Qhf GC > 0.3 enter the cosmologically interesting range 0.1 < 



< 0.3, due 



to NCC (squares with a dot in the figure), SLC (filled squares), or both NCC and SLC (empty 
squares). There are also models which enter the less interesting region for LSP to be CDM 
(flh COA < 0.1). The largest amount of models was shifted out due to NCC (filled circles), and a 
relatively small amount of models was shifted out due to SLC (circles with dots), or both NCC 
and SLC (open circles). Contrary to mSUGRA §, in the effMSSM with SLC and NCC we can 
not find a possibility to derive any cosmological upper limit for m x . There are cosmologically 
interesting LSPs within the full mass range 12 GeV < m x < 720 GeV (Fig. |J) accessible in our 
scan irrelevantly to neglecting or inclusion of any coannihilation channels in question. 
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FIG. 6. Variation of relic density against gaugino purity as function of coannihilation channels included. 
The NCC reduces the RD especially for models with higgsino-like LSPs and shifts these models out of 
cosmological interest. The NCC together with the SLC leaves only LSPs with V > 0.6 in the cosmologically 
interesting region. 

Cosmologically interesting LSPs occur with arbitrary compositions when coannihilations are 
ignored (Fig. ||), the inclusion of NCC rules out all the models with higgsino-like LSPs, and SLC 
further tends to rule out LSPs with mixed composition, so that only LSPs with V > 0.6 remain 
as dominant CDM candidates. While NCC is important both for higgsino-like and gaugino-like 
LSPs the SLC mainly affects only gaugino-like LSPs. In general our estimations (Fig. 
accordance with |8| when bino-like LSPs and SLC are concerned and in accordance with 
NCC effect is concerned. 



are m 
lOH when 



B. Detection rates 



Now we briefly consider the influence of NCC and SLC on prospects for indirect and direct 
CDM neutralino detection. Figure displays the expected indirect detection rates for upgoing 
muons produced in the Earth by neutrinos from decay products of XX annihilation which takes 
place in the core of the Earth or of the Sun. We compare the rate predictions for cosmologically 
interesting LSPs when the RD is evaluated with or without coannihilations taken into account. 
We have seen before that the RD in most models with m x < 250 GeV is untouched both by SLC 
and NCC, because the difference ttinlsp — w x is too large to yield significant effects, therefore the 
corresponding detection rates are not influenced (depicted in the figures as filled squares with dots 
inside). For XX annihilation in the Earth upgoing muon detection rates merely lie within the range 
10- 19 m~ 2 • yr" 1 < < 5 • lO^m" 2 • yr" 1 as long as m x < 250 GeV. When m x > 250 GeV, 
some of the models with 0.1 < flfif GC < 0.3 are ruled out from the cosmological interesting range 
(flh COA < 0.1; Fig. |) mainly due to NCC (Fig. [?]). Others models with f2/i 2 GC > 0.3 are shifted 
inside this region (Fig. || and Fig. 0) mainly due to SLC. In total, for m x > 250 GeV one finds 
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10 _19 m~ 2 ■ yr -1 < P^bcc < 5 • 10 _7 m~ 2 • yr -1 , when the RD is evaluated with coannihilations 
are taken into account and, 10 _19 m~ 2 • yr _1 < T^iqc < 4 • 10~ 6 m~ 2 • yr _1 when coannihilations 
are neglected. The large values of the detection rates of XX annihilation in the Earth are slightly 
decreased (from 10 _6 m -2 • yr _1 to 10~ 7 m~ 2 • yr _1 ) only for heavy LSPs m x > 500 GeV in 
accordance with the fact that the corresponding models are ruled out from the cosmologically 
interesting range. The few models with maximal detection rates at a level of 10 -4 m~ 2 • yr _1 
(m x < 500 GeV) are found to be untouched. 
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FIG. 7. Indirect detection rate for upgoing muons from \X annihilation in the Earth (a) and 
the Sun (b). As in Fig. || empty circles, grey circles, black dots, and squares correspond to 

, Q.h% cc < 0.3, respectively. NCC slightly decreases the detection rates for 



0.1 < Vth 



IGC , Q,h sljC , 0,h NCC 
models with m x > 500 GeV. 
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FIG. 8. Neutralino-proton scattering cross sections for scalar (spin-independent) interaction (a) and 
axial (spin-dependent) interaction (b). The notations as in Fig. [?]. 
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In the case of indirect detection of upgoing muons from XX annihilation in the Sun one has in 
general a similar behavior for models with m x > 250 — 300 GeV. The only noticeable difference 
is in the absolute predictions for detection rates for models with m x > 600 GeV where instead 
of r^icc < 10~ 4 m~ 2 • yr _1 one expects the rates to be T^bcc < 3 • 10~ 5 m~ 2 • yr _1 . The 
highest predicted detection rates of 10 _1 m~ 2 • yr _1 are again correlated to a few models which are 
untouched by coannihilation. 

Figure || shows neutralino-proton scattering cross sections for the scalar (spin-independent) and 
the axial (spin-dependent) interactions. As in the previous figures the models with m x < 250 GeV 
are hardly affected by coannihilation, and for the majority of those models both neutralino-proton 
and neutralino- neutron scattering cross sections reach values a < 10~ 17 GeV~ 2 with the maximal 
cross section of order 10 _15 GeV~ 2 . Cosmologically interesting models with m x > 250 GeV were 
influenced by coannihilations as discussed above, and the maximal value of the neutralino-nucleon 
cross-section decreases from 10 -15 GeV -2 to 5 • 10~ 16 GeV~ 2 for the models with m x > 500 GeV. 
In total, independently of neglection or inclusion of NCC and SLC the maximal scalar scattering 
neutralino-nucleon cross section reaches 10 _16 -10 -15 GeV -2 . 

The spin-dependent neutralino-nucleon cross sections are typically higher than the spin- 
independent ones, and we have found the maximal values at 10 -10 GeV -2 for the axial neutralino- 
proton and 10~ n GeV~ 2 for the axial neutralino-neutron scattering for the models which are 
untouched by the coannihilations. The majority of cosmologically interesting models yields axial 
neutralino-proton scattering cross sections in the range 5 • 10~ 16 GcV~ 2 < a < 2 • 10~ 12 GeV -2 
and axial neutralino-neutron scattering cross sections in the range 2 ■ 10~ 16 GeV~ 2 < a < 
8 ■ 10" 13 GeV -2 . 




100 200 300 400 500 600 700 8C 

mLSP, GeV 



FIG. 9. Event rate for direct neutralino detection in a Ge detector. As in Fig. empty circles, grey 
circles, black dots, and squares correspond to 0.1 < f2/ij GC , fi/i| LC , flh^ cc , Q.h% cc < 0.3, respectively. 
NCC slightly decreases the maximal event rates for models with m x > 500 GeV, but the models with 
smaller LSP mass are untouched by the coannihilations. 

Figure || shows the expected direct detection event rates calculated for a 73 Ge detector when 
NSS, SLC, and BCC are taken into account. For models with m x < 250 GeV coannihilations of 
any kind play no role. The optimistic estimations of the event rate for models with m x > 400 GeV 
are slightly decreased due to NCC. 
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IV. CONCLUSIONS AND OUTLOOK 



Within the low-energy effective MSSM we calculated the neutralino relic density (RD) taking 
into account both slcpton-neutralino (SLC) and neutralino-chargino/ncutralino (NCC) coannihi- 
lation channels. To this end we have implemented in our code JlO"[ the relic density part (with 
neutralino-chargino coannihilations) of the DarkSusy code jyj supplied with the adopted code of 
H (calculating slepton- neutralino coannohilations). 

We have shown that in effMSSM the maximum factors of RD decrease due to NCC as well 
as due to SLC can reach 10~ 3 , as long as the lower limits for rrif and m^± are similar (both 
being of order of 80 GeV). We conclude that both coannigilation effects are comparable in the 
effMSSM. For the majority of models affected by coannihilations and successfully passed all relevant 
accelerator, cosmological and rare-decay constraints it was observed that either NCC or SLC alone 
produces significant reduction of RD while the other coannihilation channel gives considerably 
smaller or zero reduction effect. Both coannihilations (NCC and/or SLC) are found to produce 
non-negligible effect only if the relevant NLSP mass is smaller than 1.15m x . The type of the 
NLSP determines the dominant coannihilation channel. Due to the fact that the effMSSM more 
often favors neutralino and chargino, but not sleptons to be the NLSP (the NLSP-LSP mass 
differences in general are systematically larger for sleptons than for gauginos) the NCC is the more 
often dominant coannihilation channel in agreement with JTcj ] . Only LSPs with purity V > 0.6 
remain CDM candidates of cosmological interest. Some models with Q,h 2 > 0.3 under neglection of 
coannihilation enter into the cosmologically interesting region merely due to SLC, and some other 
models shift out of the region below flh 2 < 0.1, merely due to NCC. In the effMSSM, contrary to 
mSUGRA H , both coannihilations do not imply new cosmological limits on the mass of the LSP. 
We noticed that the most optimistic predictions for neutralino-nucleon cross sections, indirect and 
direct detection rates for cosmologically interesting models are untouched by these coannihilations. 
Only for large m x > 500 GeV, the respectively high values are slightly reduced, because of the 
NCC rules out corresponding models from the cosmological interesting region 0.1 < flh 2 GC < 0.3. 

When our paper was almost finished we saw the preprints [^| J3^J3^ |. In a new sophisti- 
cated code micrOMEGAs for calculations of the relic density in the MSSM is presented. The main 
characteristics of this code includes complete tree-level matrix elements for all subprocesses; all 
coannihilation channels with neutralinos, charginos, sleptons, squarks and gluinos; loop-corrected 
Higgs masses and widths. All calculations are performed with CompHEP fljsfl . In |37j the relic den- 
sity of neutralinos in the mSUGRA was calculated on the basis of annihilation diagrams, involving 
sleptons, charginos, neutralinos and third generation squarks. The code CompHEP Q was used, 
too, here. It was found in |37]] that coannihilation effects are only important on the edges of the 
model parameter space, where some amount of fine-tuning is necessary to obtain a reasonable relic 
density. This paper is mostly aimed at prospects for SUSY search with various e + e~ and hadron 
colliders and pays no attention on the interplay between different coannihilation channels. 

In addition, a full set of exact, analytic expressions for the annihilation of the lightest neutralino 
pairs into all two-body tree-level final states in the framework of minimal SUSY is now available 
]3q] . This set of expressions does not rely on the partial wave expansion, includes all the terms 
and is valid both near and further away from resonances and thresholds for new final states. 
Further extension of this approach to coannihilation processes together with the above-mentioned 
CompHEP-based codes will supply perhaps one of most powerful tools for complete relic density 
calculations. 

V.B. thanks the Max Planck Institut fuer Kernphysik for the hospitality and RFBR (Grant 
00-02-17587) for support. 
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